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■ Abstract. Echo mapping exploits light travel time delays, revealed by 

Q , multi- wavelength variability studies, to map the geometry, kinematics, 

and physical conditions of reprocessing sites in photo-ionized gas flows. 
In AGNs, the ultra-violet to near infra-red light arises in part from repro- 
cessing of EUV and X-ray light from a compact and erratically variable 
source in the nucleus. The observed time delays, 0.1-2 days for the con- 
tinuum, 1-100 days for the broad emission lines, probe regions only micro- 
arcseconds away from the nucleus. The continuum time delays map the 
temperature-radius profiles of the AGN accretion discs. The emission- 
£vq ■ line delays reveal radially stratified ionization zones, identify the nature 

of the gas motions, and estimate the masses of the central black holes. 
By using light travel time to measure the sizes of AGN accretion discs 
and photo-ionized zones, echo mapping offers two independent ways to 
^ ■ measure redshift-independent distances to AGNs. 

6 

52 ' 1. Introduction 

C3 ■ 

Following this brief introduction, section |2] reviews in general terms the tech- 
k> | niques used for echo mapping of active galaxies. Sections ^J, [4], and |5] then 

5_i ■ discuss the use of echo mapping experiments to probe the geometry, the kine- 

matics, and the physical conditions in Seyfert 1 galaxies. Section ^] considers 
two direct methods based on echo mapping to determine redshift-independent 
distances, and hence cosmological parameters Hq and qo. 

1.1. The Black Hole Accretion Disc Model for AGNs 

The standard model of an Active Galactic Nucleus (AGN) envisions a super- 
massive (10 6_9 Mq) black hole in the core of a galaxy. The black hole is fed gas 
from an accretion disc. Around the disc is a geometrically thick torus of dust 
and molecular gas. The disc and torus receive gas from the galaxy's interstellar 
medium, and from various processes (tides, bow shocks, winds, irradiation) that 
strip material from the stars that pass through this region. 

The ultra-violet and optical spectra of quasars, and their lower-luminosity 
cousins in the nuclei Seyfert galaxies, have two types of emission lines: narrow 
forbidden and permitted emission lines (v £ 1000 km s _1 ), and broad permitted 
emission lines (v ^ 10,000 km s^ 1 ). The narrow lines are constant; the contin- 
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uum and the broad lines vary. The narrow lines arise from lower-density gas at 
larger distances from the nucleus. 

Seyfert nuclei come in two types thought to represent different viewing 
angles. Seyfert 1 spectra exhibit both broad and narrow emission lines, while 
in Seyfert 2 spectra the broad lines are absent or visible only in a faint linearly 
polarized component of the spectrum. The current consensus is that the thick 
dusty torus blocks our view of the BLR when viewed at i ^ 60° , but some BLR 
light scatters toward us after rising far enough above the plane to clear the torus. 

The Narrow Line Region (NLR) in nearby Seyfert galaxies is resolved by 
HST imaging studies with resolutions of ~ 0.1 arcseconds. The NLR typically 
has a clumpy bi-conical morphology, aligned with the axis of radio jets, and sug- 
gesting broadly-collimated ionization cones emerging from the nucleus, perhaps 
collimated by the dusty torus. The BLR and continuum production regions are 
unresolved by HST. Echo mapping experiments use light travel time delays re- 
vealed by variability on 1-100 day timescales to probe these regions on ~ 10~ 6 
arcsecond scales. 

1.2. Photo-Ionization Models 

Just how large is the BLR? If the emission lines are powered by photo-ionization, 
then we may employ a model to estimate the size of the photo-ionized region. 
Photo-ionization models such as CLOUDY (Ferland, et al. 1998) consider the 
energy and ionization balance inside a 1-dimensional gas cloud parameterized 
by hydrogen number density bh, column density A^h, and distance R from the 
source of ionizing radiation L\. The calculated emission- line spectrum emerging 
from such a gas cloud depends primarily on the ionization parameter, 

ionizing photons Q 

target atoms 47i\R 2 cnH' 

where the luminosity of hydrogen-ionizing photons is 

with L\ the luminosity spectrum emitted by the nucleus. Re-arranging this 
equation yields 

1/2 

R 



(-Q—Y (3) 

\A7rUcnuJ ' 



By comparing the flux ratios observed for the broad emission lines in the spectra 
of AGNs against those predicted by the single-cloud photo-ionization models, 
typical parameters in the photo-ionized gas are found to be U ~ 10~ 2 and 
ran ~ 10 8_10 cm~ 3 . With these conditions, the light travel time across the radius 
of the ionized zone is 

R I Q \V 2 / U V V V "H V' /2 



. 20M llP^J llF^J 1 10^ ) • (4) 

where Q ~ 10 54 s _1 (iJo/100)~ 2 is appropriate for the Seyfert 1 nucleus of NGC 5548. 
On the basis of such predictions, astronomers searched for reverberation effects 
on timescales of months before eventually realizing that reverberation was oc- 
curring on much shorter timescales. 
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Figure 1. Ionizing photons from a compact source are reprocessed 
by gas clouds in a thin spherical shell. A distant observer sees the re- 
processed photons arrive with a time delay ranging from 0, for the near 
edge of the shell, to 2R/c, for the far edge. The iso-delay paraboloids 
slice the shell into zones with areas proportional to the range of delays. 



2. Echo Mapping Methods 
2.1. Reverberation 

A compact variable source of ionizing radiation launches spherical waves of heat- 
ing or cooling that expand at the speed of light, triggering changes in the radi- 
ation emitted by surrounding gas clouds. The light travel time from nucleus to 
reprocessing site to observer is longer than that for the direct path from nucleus 
to observer. Thus a distant observer sees the reprocessed radiation arrive with 
a time delay reflecting its position within the source. The time delay for a gas 
cloud located at a distance R from the nucleus is 

R 

r= _(l + cos0) , (5) 
c 

where the angle 9 is measured from for a cloud on the far side of the nucleus 
to 180° on the line of sight between the nucleus and the observer. The delay is 
for gas on the line of sight between us and the nucleus, and 2R/c for gas directly 
behind the nucleus (Fig. |]). The iso-delay contours are concentric paraboloids 
wrapped around the line of sight. Such time delays are the basis of echo mapping 
techniques. 

The validity of this basic picture is well supported by the results of inten- 
sive campaigns designed to monitor the variable optical, ultra-violet, and X-ray 
spectra of Seyfert 1 galaxies. Notable among these are the AGN Watch cam- 
paigns (http://www.astronomy.ohio-state.edu/~agnwatch/). In Seyfert 1 
galaxies, the variations observed in continuum light are practically simultaneous 
at all wavelengths throughout the ultra-violet and optical, but the correspond- 
ing emission line variations exhibit time delays of 1 to 100 days, probing the 
photo-ionized gas 1 to 100 light days from the nucleus. This corresponds to 
micro- arcsecond scales in nearby Seyfert galaxies. 



3 



2.2. Linear and Non-Linear Reprocessing Models 

The ionizing radiation that drives the line emission includes X-ray and EUV light 
that cannot be directly observed. However, quite similar variations are seen in 
continuum light curves at ultra-violet and optical wavelengths. These observable 
continuum light curves provide suitable surrogates for the unobservable ionizing 
radiation. The reprocessing time (hours) is small compared with light travel 
time (days) in AGNs. 

In the simplest reprocessing model, the line emission L(t) is taken to be 
proportional to the continuum emission C(t — t), with a time delay r arising 
from light travel time and from local reprocessing time, thus 

L(t) = *C(t - r) . (6) 

Since the reprocessing sites span a range of time delays, the line light curve L(t) 
is a sum of many time-delayed copies of the continuum light curve C (t) , 

poo 

L(t) = / *(t) C{t-r) dr. (7) 
Jo 

This convolution integral introduces VP(t), the "transfer function" or "convolu- 
tion kernel" or "delay map" , to describe the strength of the reprocessed light 
that arrives with various time delays. 

Since light travels at a constant speed, the surfaces of constant time delay 
are ellipsoids with one focus at the nucleus and the other at the observer. Near 
the nucleus, these ellipsoids are effectively paraboloids (Fig. ||). ^(r) is in effect 
a 1-dimensional map that slices up the emission-line gas, revealing how much 
gas is present between each of the iso-delay paraboloids. 

The aim of echo mapping is to recover ^(t) from measurements of L(t) and 
C(t) made at specific times ij. To fit such observations, the linear model above 
is too simple in at least two respects. The first problem is additional sources of 
light contributing to the observed continuum and emission- line fluxes. Examples 
are background starlight, and narrow emission lines. When these sources do not 
vary on human timescales, they simply add constants to L(t) and C(t), 

POO 

L(t) = L+ #(r) [C{t -t)-C] dr. (8) 
Jo 

A second problem with the linear reprocessing model is that the reprocessed 
emission is more generally a non-linear function of the ionizing radiation. For 
example, when a cloud is only partially ionized, the line emission increases with 
the ionizing radiation, but as the cloud becomes fully ionized the line emission 
may saturate or even decrease with further increases in ionizing radiation. For 
single clouds we therefore expect the response function L{C) to be highly non- 
linear (Fig. [2]). The marginal response dL/dC is generally less than the mean 
response L/C, and may become negative when an increase in ionizing radiation 
reduces the line emission. 

The total response is of course a sum of responses from many different gas 
clouds, those closer to the nucleus being more fully ionized than those farther 
away. If clouds are present at many radii, so that zones of constant ionization 



4 



Non-Linear Response of Single Cloud 



Sum of Many Non — Linear Clouds 




Figure 2. The emission-line response from gas clouds exposed to dif- 
ferent levels of photo-ionizing continuum radiation. Left: The initially 
linear response saturates and may even decrease as the cloud becomes 
more fully ionized. Right: When gas clouds are present at many radii, 
the ionization zone can simply expand, resulting in an ensemble re- 
sponse that increases monotonically and is less non-linear than that 
from a single cloud. The tangent line is then a useful approximation. 



parameter can simply expand, then the effect of averaging over many clouds is 
to produce a total response that is monotonically increasing and less strongly 
non-linear. Such a response may be adequately approximated by a tangent line, 

L(C) = L + ^(C-C) , (9) 

for ionizing radiation changes in some range above and below a mean level 
(Fig. ||). Observations support the use of this approximation - plots of ob- 
served line fluxes against continuum fluxes show roughly linear relationships, 
with L/C > dL/dC so that extrapolation to zero continuum flux leaves a posi- 
tive residual line flux. 

For these reasons, it is usually appropriate in echo mapping to adopt the 
linearized reprocessing model, equation ||, with the "background" fluxes, L for 
the line and C for the continuum, set somewhere in the range of values spanned 
by the observations. In this model, the delay map ^(t) senses each gas cloud 
in proportion to its marginal response to a change in ionizing radiation. The 
roughly linear responses from partially-ionized clouds are fully registered, while 
the saturated or diminished responses of more fully ionized clouds have a reduced 
effect. 

2.3. Cross-Correlation Analyses 

With reasonably complete light curves, it is usually obvious that the highs and 
lows in the line light curves occur later than those in the continuum. To quantify 
this time delay, or lag, a common practice is to cross-correlate the line and 
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NGC 7469 Light Curves Cross-Correlation Functions 




J I I I I I I I I I ._! I I I I I I I I I I I I L 



250 275 -10 10 

Julian date (-2450000) Time delay (days) 

Figure 3. The left-hand columns show light curves of NGC 7469 
obtained with IUE during an intensive AGN Watch monitoring cam- 
paign during the summer of 1996. The right-hand column shows the 
result of cross-correlating the light curve immediately to the left with 
the 1315 A light curve at the top of the left column; the panel at the top 
of the right column thus shows the 1315 A continuum auto-correlation 
function. Data from Wanders, et al. (1997). 
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continuum light curves. The cross-correlation function (CCF) is 



L*C(t) = / L{t) C(t-r) dt 



(10) 



Several methods have been developed and refined to compute CCFs from noisy 
measurements available only at discrete unevenly-spaced times, either by inter- 
polating the data (Gaskell & Sparke 1986, White & Peterson 1994) or binning 
the CCF (Edelson & Krolik 1988, Alexander 1997). The resulting CCFs gener- 
ally have a peak shifted away from zero in a direction indicating that changes 
in the emission lines lag behind those in the continuum (e.g. Fig. |J). The CCF 
lag - tccf _ at which L ★ C has its maximum value, serves to quantify roughly 
the size of the emission-line region. 

Since a range of time delays is present, described by ^f(r), what delay 
is measured by the cross-correlation peak? Since L(t) is itself a convolution 
between C(t) and \£(t), and since convolution is a linear operation, we have 



The CCF is therefore a convolution of the delay map with the continuum auto- 
correlation function (ACF), 



The ACF is symmetric in r, and thus always has a peak at r = 0. With rapid 
continuum variations, the ACF is sharp, and the CCF peak should be close to 
the strongest peak of \I/(t). This tends to favour short delays from the inner 
regions of the BLR (Perez, et al. 1992a). When continuum variations are slow, 
however, the ACF is broad, smearing out sharp peaks in vE'(t), and shifting 
tccf toward the centroid of ^(t). Thus the cross-correlation peak depends not 
only on the delay structure in vI/(t), but also on the character of the continuum 
variations. Different observing campaigns may therefore yield different lags even 
when the underlying delay map is the same. 

2.4. Three Echo Mapping Methods 

More refined echo mapping analyses aim to recover the delay map ^(r) rather 
than just a characteristic time lag. These echo mapping methods generally 
require more complete data than the cross-correlation analyses. Three practical 
methods have been developed. 

The Regularized Linear Inversion method (RLI) (Vio, et al. 1994, Kro- 
lik & Done 1995) notes that the convolution integral, 



L * C = (* * C) * C = * * (C * C). 



(11) 




(12) 




(13) 



becomes a matrix equation, 



L{t i )=Y J C{t i -T 3 ) *(Tj) dr 



(14) 



3 



7 



when the times tj are evenly spaced. If the times are not evenly spaced, interpo- 
late. If the matrix Qj = C(ti — tj) can be inverted, solving the matrix equation 
yields 

nrk) = Y, C ki L {ti)/dr. (15) 
k 

If the continuum variations are unsuitable, the matrix has small eigenvalues and 
the inversion is unstable, strongly amplifying noise in the measurements L(ti) 
and C(ij). This problem is treated by altering the matrix to reduce the influence 
of small eigenvalues, thereby reducing noise but blurring the delay map. 

The related method of Subtractive Optimally-Localized Averages 
(SOLA) (Pijpers & Wanders 1994) aims to estimate the delay map as weighted 
averages of the emission line measurements, 

$( r ) = J K(t, t) L(t) dt. (16) 

Since L = C * we may write this as 

$(r) = J K(t, t) J C(t - s) #(s) ds dt . (17) 

The estimate * is therefore a blurred version of the true delay map With 
suitable continuum variations, the weights K(t, t) can be chosen to make the 
blur kernel 

b( T , s) = J K(t, t) C(t - s) dt (18) 

resemble a narrow Gaussian of width A centred at s = r. The parameter A 
then controls the trade-off between noise and resolution in reconstructing *(r). 
These direct inversion methods work best when the observed light curves detect 
continuum variations on a wide range of timescales. 

The Maximum Entropy Method (MEM) (Home, Welsh & Peterson 
1991, Home 1994) is a very general fitting method allowing the use of any 
linear or non-linear reverberation model. As an extension of maximum likelihood 
techniques, MEM employs a "badness of fit" statistic, x 2 > to judge whether 
the echo model being considered achieves a satisfactory fit to the data. The 
requirement % 2 /N ~ 1 ± y/2/N, where N is the number of continuum and line 
measurements, ensures that model fits as well as is warranted by the error bars 
on the data points, without over-fitting to noise. MEM fits are required also 
to maximize the 'entropy', which is designed to measure the 'simplicity' of the 
map. For positive maps pi > 0, the entropy 

S = ^pi- qt-pilnpi/qi (19) 

i 

is maximized when pi = qi. MEM thus steers map pi toward default values 
With the default map set to 

Qi = (Pi-iPi+i) 1/2 , (20) 
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Figure 4. Simulation test of the maximum entropy method showing 
the recovery of a delay map (top left) from data points sampling an er- 
ratically varying continuum light curve (bottom) and the corresponding 
delay-smeared emission- line light curve (top right). The reconstructed 
delay map closely resembles the true map ^(t) oc re~ T . 

the entropy steers each pixel toward its neighbors, and MEM then finds the 
'smoothest' positive map that fits the data. For further technical details see 
Home (1994). 

Fig. ^| illustrates recovery of a delay map from a MEM fit to fake light 
curves. The lower panel shows an erratically varying continuum light curve. 
The line light curve in the upper panel is smoother and has time-delayed peaks. 
A smooth continuum light curve C{t) threads through the data points, and 
extrapolates to earlier times. Convolving this light curve with the delay map 
^(t), shown in the upper left panel, gives the line light curve L(t), fitting the 
data points in the upper right panel. Dashed lines give the backgrounds L and 
C. The MEM fit with x 2 /N = 1 adjusts C{t), ^(r), and L to fit the data points 
while keeping C(t) and ^(t) as smooth as possible. 

In this simulation test, the true transfer function, ^(t) oc re~ r , is accurately 
recovered from the data. The fake dataset represents the type of data that could 
be obtained with daily sampling over a baseline of 1 year. While this is rather 
better than has been achieved so far, future experiments (Kronos, Robonet) 
specifically designed for long-term monitoring will make this simulation more 
relevant. 

3. Mapping the Geometry of Emission-Line Regions 
3.1. Spherical Shells 

Having shown that we can recover the delay map ^(t), given suitable data, 
how may we interpret it? This is relatively straightforward if the geometry is 
spherically symmetric. 
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Figure 5. Delay maps for spherical geometry. Left: When the 
line emission is isotropic, the contribution from each spherical shell 
is constant from to 2R/c. The total delay map, summing many 
spherical shells, must decrease monotonically. Right: When the line 
emission is directed inward, toward the nucleus, the response at small 
delays is reduced, so that each shell's contribution increases with r 
between and 2R/c. The total delay map then rises to a peak away 
from zero. 



Consider a thin spherical shell that is irradated by a brief flash of ionizing 
radiation from a source located at the shell's centre. The flash reaches every 
point on the shell after a time R/c. Since recombination times are short, each 
point responds by emitting a brief flash of recombination radiation. A distant 
observer sees first the flash of ionizing radiation, and then the response of repro- 
cessed light from the shell arriving with a range of time delays. The time delay 
is for the near edge of the shell, and 2R/c on the far edge of the shell. The 
iso-delay paraboloids slice up the spherical shell into zones with areas propor- 
tional to the range of time delay (Fig. [I]). The response from the spherical shell 
is therefore a boxcar function, constant between the delay limits and 2R/c 
(Fig. §. 

Any spherically symmetric geometry is just a nested set of concentric spher- 
ical shells. The delay map for a spherical geometry is therefore a sum of boxcar 
functions (Fig. ||). Since all the boxcars begin at delay 0, the delay map must 
peak at delay 0, and decrease monotonically thereafter. The contribution from 
the shell of radius R may be identified from the slope d^/dr evaluated at the 
appropriat e tim e delay r = 2R/c. (The case of anisotropic emission is discussed 



in Section 3.3. below. 



3.2. The BLR is Smaller than Expected 

Fig. [6| shows a maximum entropy fit of the linearized echo model of Eqn. || to 
H/3 and optical continuum light curves of NGC 5548 (Home, Welsh & Peterson 
1991). Subtracting the continuum background C, convolving with the delay 
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Figure 6. Echo maps of H/3 emission in NGC 5548 found by a max- 
imum entropy fit to data points from a 9-month AGN Watch monitor- 
ing campaign during 1989. The optical continuum light curve (lower 
panel) is convolved with the delay map (top left) to produce the H/3 
emission line light curve (top right). Horizontal dashed lines give the 
mean line and continuum fluxes. Three fits are shown to indicate likely 
uncertainties. Data from Horne, Welsh & Peterson (1991). 



map ^(t), and adding the line background L, gives the H/3 light curve. The 
three fits shown, all with x 2 /N = 1, indicate the likely range of uncertainty due 
to the trade-off between the 'stiffness' of ^(r) and C{t). 

The H/3 map has a single peak at a delay of 20 days, and declines to low 
values by 40 days. This suggests that the size of the H/3 emission-line region 
is 10-20 light days. This is 10 to 20 times smaller than the 200 light day size 
estimated in Eqn. ||| on the basis of single-cloud photo-ionization models. With 
clouds this close to the nucleus, the ionization parameter U would be higher 
unless gas densities are increased by a factor of 100, to kh ~ 10 n cm~ 3 . 

3.3. Anisotropic Emission 

The H/3 response in NGC 5548 is smaller at delay than at 20 days. This lack 
of a prompt response conflicts with the monotonically decreasing delay map we 
expect for a spherical geometry. One interpretation is that there is a deficit of 
gas on the line of sight to the nucleus. 

However, more likely this is a signature of anisotropic emission of the H/3 
photons arising from optically thick gas clouds that are photo-ionized only on 
their inward faces (Ferland, et al. 1992, O'Brien, et al. 1994). The preference 
for inward emission of H/3 photons reduces the prompt response by making the 
reprocessed light from clouds on the far side stronger than that from clouds on 
the near side of the nucleus. With inward anisotropy, the impulse response of a 
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spherical shell is a wedge, increasing with r, rather than a boxcar (Fig. |5|). The 
sum of wedges has a peak away from zero. 

3.4. Stratified Temperature and Ionization 

Ultra-violet spectra from IUE and HST provide shorter-wavelength continuum 
light curves, and light curves for a variety of high and low ionization emission 
lines. These are useful probes of the temperature and ionization structure of the 
reprocessing gas. 

The continuum light curves display very similar structure at all optical and 
ultraviolet wavelengths. The continuum lightcurve at the shortest ultra-violet 
wavelength observed is normally used as the light curve against which to study 
time delays at other continuum wavelengths and in the emission lines. In most 
cases studied to date the optical and ultra-violet continuum light curves exhibit 
practically simultaneous variations (e.g. Peterson, et al. 1998), implying that 
the continuum production regions are unresolved, smaller than a few light days. 

In one case, NGC 7469, delays of 0.1 to 1.7 days are detected between 
the 1350A and 7500A continuum light curves based on 40 days of continuous 
IUE and optical monitoring (Wanders et al. 1997, Collier et al. 1998). The 
delay increases with wavelength as r a A 4 / 3 , suggesting that the temperature 
decreases as T oc i?~ 3//4 . 

A systematic pattern is generally seen in the emission- line time delays. 
High-ionization lines (N v, He n) exhibit the shortest delays, while lower-ionization 
lines (C iV,Lya,H/3) have longer delays. This was established in the first AGN Watch 
campaign, by cross-correlation analysis (Clavel et al. 1991) and by maximum 
entropy fits (Krolik, et al. 1991), using IUE and optical light curves of NGC 5548 
sampled at 4 day intervals for 240 days. The effect is now seen in many objects 
(Table ||). This pattern is consistent with photo-ionization models, in which 
higher ionization zones occur closer to the nucleus. 



Table 1. AGN Watch Monitoring Campaigns 

target year tccf (days) references 

Hen Lya Civ H/3 



NGC 5548 


1989 


7 


12 


12 


20 


1, 2, 3 


NGC 3783 


1992 





4 


4 


8 


4, 5 


NGC 5548 


1993 


2 


8 


8 


14 


6 


Fairall 9 


1994 


4 


17 




23 


7,8 


3C 390.3 


1995 


10 


50 


50 


20 


11, 12 


NGC 7469 


1996 


1 


3 


3 


6 


9, 10 



References: 1. Clavel, et al. 1991; 2. Peterson, et al. 1991; 3. Dietrich, 
et al. 1993; 6. Korista, et al. 1995; 4. Reichert, et al. 1994; 5. Stirpe, et al. 1994; 
7. Rodriguez-Pascual, et al. 1997; 8. Santos-Lleo, et al. 1997; 9. O'Brien, 
et al. 1998; 10. Dietrich, et al. 1998; 11. Wanders, et al. 1997; 12. Collier, 
et al. 1998. 

Another interesting result is that the time delay for Civ emission is gen- 
erally smaller than that for C III] emission. In early work, this line ratio was 
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Figure 7. The radii of ionized zones emitting broad H/3 emis- 
sion lines for AGNs with different ionizing continuum luminosities. 
The radii are estimated from time delays that are found by cross- 
correlating H/3 and optical continuum light curves. The diagonal line 
with R oc L 1 / 2 corresponds to a constant ionization parameter. Figure 
from Kaspi, et al. (1996). 

the basis for estimating the gas density in the BLR to be bh ~ 

10 8-io 

cm . 

However, this practice is now considered dubious because their different time 
delays indicate that these lines arise in different regions, The Civ line arises 
from gas clouds closer to the nucleus, and a higher density, bh ~ 10 11 cm -3 , is 
therefore needed to maintain the ionization parameter. 

3.5. Ionized Zones Expanding with Luminosity 

The ionization zones should be larger in higher luminosity objects. This predic- 
tion has been tested using echo mapping studies of a dozen active galaxies, in- 
cluding two low-luminosity quasars, spanning three decades in luminosity (Fig. |7[ 
Kaspi, et al. 1996). The lag for each object is determined by cross-correlating 
H/3 and optical continuum light curves. A log-log plot of the time delay vs lu- 
minosity shows a standard deviation of about 0.3 dex about a correlation of the 
form 

th/3 ~ 17d L\{ 2 , (21) 

where L44 is the hydrogen-ionizing luminosity in units of 10 44 erg s _1 . For com- 
parison, if gas clouds with similar densities are present over a wide range of ra- 
dius, then zones of constant ionization parameter should increase with R oc L 1 / 2 . 

In the above correlation, derived by comparing different objects, the ob- 
served scatter may be due in part to differences in inclination, black hole mass, 
and accretion rate. The ionized zone in each object is expected expected to 
change size if the ionizing luminosity changes by a substantial factor. There is 
some evidence for this, e.g. from 8 years of monitoring NGC 5548 (Peterson, 
et al. 1999). To map this effect, the echo model may be generalized to allow the 
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Figure 8. Velocity-delay maps for spherical freefall and Keplerian 
disc kinematics. Spherical shells map to diagonal lines, and disc annuli 
map to ellipses. The central mass is 1O 7 M in both cases. The disc 
inclination angle is 60°. 

delay map to change with the driving continuum light curve, 



4. Mapping the Kinematics of Emission-Line Regions 

Combining time delays from variability with Doppler shifts from emission-line 
profiles yields a velocity-delay map ty(v,r) to probe the kinematics of the flow. 
Fig. U illustrates this capability by sketching the formation of velocity-delay maps 
for spherical free-fall into a point-mass potential, and for a Keplerian disc (Welsh 
& Home 1991, Perez, et al. 1992b). The dramatically different appearance of the 
two flows highlights the power of velocity-delay maps for kinematic diagnosis. 

A radial flow, inward or outward, produces a strong red-blue asymmetry in 
the velocity-delay map, while an azimuthal flow does not. Line emission from 
gas flowing in toward the nucleus is redshifted (v > 0) on the near side (r small) 
and blueshifted (v < 0) on the far side (r large). The signature of inflowing 
gas is therefore small time delays on the red side and large time delays on the 
blue side. Outflowing gas produces just the opposite red-blue asymmetry. Gas 
circulating around the nucleus has the same time delay on the red and blue side, 
producing a symmetric velocity-delay map. 

Any spherically-symmetric flow is a nested set of thin spherical shells. The 
time delay r and Doppler shift v are 




(22) 



r = 



f (1 + COS0) 



v = vr cos 9 



(23) 
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with R the shell radius, 9 the angle from the back edge, and vr the outflow 
velocity. The linear dependence of both t and v on cos 9 implies that each shell 
maps into a diagonal line in the velocity- delay plane, as shown in Fig. ||. 

The circulating disc flow is a set of concentric co-planar cylindrical annuli. 
The time delay and Doppler shift are 

t = — (1 + sinisin^) , v = s'mi cos 4>, (24) 

with R and i the radius and inclination of the annulus, <f> the azimuth, and iu 
the azimuthal velocity. Each annulus maps to an ellipse on the velocity-delay 
plane, as in Fig. |8[ Inner annuli map to "squashed" ellipses (large ±t>, small r), 
while outer annuli map to "stretched" ellipses (small ±u, large t). 

While $>(v,t) is a distorted 2-dimensional projection of the 6-dimensional 
phase space, it can reveal important aspects of the flow, particularly if the 
velocity field is ordered and has some degree of symmetry. The envelope of 
the velocity-delay map may reveal the presence of virial motions v 2 oc GM/ct, 
with M the mass of the central object. A red/blue asymmetry, or lack thereof, 
gauges the relative importance of radial and azimuthal motions. Disordered 
velocity fields with a range of velocities at each position smear the map in the v 
direction. The far side (large r) may be enhanced by anisotropic emission from 
optically thick clouds radiating their lines inward toward the nucleus (Ferland, 
et al. 1992, O'Brien, et al. 1994). 

To derive Doppler-delay maps from the data, we simply slice the observed 
line profile into wavelength bins, and recover the delay map at each wavelength. 
This is a simple extension of the echo model used to fit continuum and emission- 
line light curves. At each wavelength A and time t the emission-line flux 

POO 

L{\,t) = L{\)+ [C(t-T)-C] #(A,r) dr (25) 
Jo 

is obtained by summing time-delayed responses to the continuum light curve 
and adding a time-independent background spectrum L(X). 

We observe C(t) and L(\,t) at some set of times U. We fix C, e.g. close 
to the mean of the observed values. We then adjust the continuum light curve 
C{t), the background spectrum L(X), and the Doppler-delay map ^(A,t), in 
order to fit the observations. The fitting procedure maximizes the entropy to 
find the "simplest" map(s) that fit the observations with x 2 /N = 1. 

Fig. |9] exhibits the velocity-delay map for Civ 1550 and Hen 1640 emission 
as reconstructed from 44 IUE spectra of NGC 4151 covering 22 epochs during 
1991 Nov 9 - Dec 15 (Ulrich &: Home 1996). While spanning only 36 days, 
this campaign recorded favorable continuum variations, including a bumpy ex- 
ponential decline followed by a rapid rise, which were sufficient to support echo 
mapping on delays from to 20 days. 

Disregard the strong C iv absorption feature, which obliterates the delay 
structure at small velocities. The wide range of velocities at small delays and 
smaller range at larger delays suggests virial motions. The dashed curves in 
Fig. H give escape velocity envelopes v = \J2GM /cr for masses 0.5, 1.0, and 
2.0 x 10 7 M . A mass of order 10 7 M Q may be concentrated within 1 light-day 
of the nucleus. 
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Figure 9. memecho velocity-delay map of Civ 1550 emission (and 
superimposed Hen 1640 emission) in NGC 4151. Dashed curves give 
escape velocity for masses 0.5, 1.0, and 2.0 x 1O 7 M . 

The approximate red-blue symmetry of the C IV map rules out purely radial 
inflow or outflow kinematics. The somewhat stronger C iv response on the red 
side at small delays and on the blue side at larger delays suggests a gas com- 
ponent with freefall kinematics. However, if the C iv emission arises from the 
irradiated faces of optically-thick clouds, those on the far side of the nucleus will 
be brighter, and the redward asymmetry can then be interpreted as an outflow 
combined with the inward C IV anisotropy. 

Velocity-delay maps have so far been constructed only for C iv emission 
in NGC 5548 (Wanders, et al. 1995, Done & Krolik 1996) and NGC 4151 
(Ulrich & Home 1996). Both systems show the same trend of velocity dispersion 
decreasing with delay, and red response stronger at small delays. As more maps 
are constructed, we will learn whether these are general characteristics of Seyfert 
broad-line regions. 

5. Mapping Physical Conditions in Emission-Line Regions 
5.1. Quasar Tomography 

The rich information coded in high quality time-resolved spectrophotometry of 
a reverberating emission-line spectrum can be extracted only by fitting observa- 
tions in far greater detail than has previously been attempted. Here we extend 
previous echo mapping methods to map simultaneously the geometry, kinemat- 
ics, and physical conditions characterizing the population of photo-ionized gas 
clouds. To do this we fit simultaneously the complete reverberating spectrum, 
including the reverberations observed in the fluxes and profiles of numerous emis- 
sion lines. These fits incorporate explicitly the predictions of a photo-ionization 
model such as CLOUDY, thus accounting for the non-linear and anisotropic 
responses of emission-line clouds to changes in the ionizing radiation. 
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We characterize a gas cloud by its density nn, column density Ah, distance 
from the nucleus R, azimuth 0, and Doppler shift v. The entire population of 
gas clouds is then described by a 5-dimensional map ^(R,9,nu, Nn,v). We 
omit the cloud's position angle <j) around the line of sight, and the perpendicular 
velocity components, because the data provide no useful constraint on these. 

Assuming that the shape of the ionizing spectrum is known, and that the 
time-dependent ionizing photon luminosity is Q(t), then the ionization param- 
eter obtaining at time t is 

U{t) = W~ T) . (26) 

The reprocessing efficiency €l{U, nn, Nr, 0) for emission line L depends in a 
unique way on U, nn, and Nr. The reprocessing efficiencies are evaluated 
with a photo-ionization code, e.g. CLOUDY, for chosen element abundances. 
To save computer time, these are pre-calculated on a suitable grid of cloud 
parameters, e.g. equally spaced in logU, lognn, and log Ah (Korista, et al. 
1997). Results required for any cloud parameters are subsequently found rapidly 
by interpolation in the grid. 

Because clouds are irradiated on their inward faces, the reprocessed line 
emission is generally anisotropic. We allow for this by letting depend on 9. 
CLOUDY computes the emission emerging inward, = 0, and outward, 9 = 90°. 
For intermediate angles we interpolate linearly in cos0. 

The observed spectrum is a sum of three components: direct light from 
the nucleus, reprocessed light from the surrounding gas clouds, and background 
light, - e.g. from stars, 

U{\ t) = /,f (A, t) + /* (A, t) + f? (A) . (27) 

The direct light is 

f?M = 9 <») 

where D is the distance and the spectral shape is 

«,(A) = MM>. (29, 

The reprocessed light is a sum over numerous emission lines, where Al and 
€l are the rest wavelength and the reprocessing efficiency for line L. The gas 
clouds are distributed over a 3-dimensional volume. At each location the cloud 
population has a distribution over density njj, column density Ah, and Doppler 
shift v. The reprocessed light arising from such a configuration is 



fv(\, t)= J 2-irRdR dn dn u dN H dv V(R, 9, n H , N u ,v) 

EL^(U,n R ,N u ,9)G u (X-(l + v/c)X L ) (30) 



Here fi = cos0, G v is a Gaussian profile to apply the appropriate Doppler shift, 
and the Dirac 5 functions ensure that the appropriate time delay and ionization 
parameter are used at each reprocessing site. 
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Figure 10. Reconstructed map 0,nn) of a thin spherical shell 
(left) from a maximum entropy fit to 7 ultra-violet emission-line light 
curves (right). The light curve fits achieve x 1 /N = 1. The shell radius 
R/c = 12d and density log ran = 11 are correctly recovered. For further 
details, see text. 



To fit the above model to observations of f v (X, t), we adjust D, Q(t), f„(\), 
and ^f(R, 0, nn, Nn,v). The 5-dimensional cloud map \P can have loads of pixels, 
~ 10 5-6 . The constraints available from reverberating emission-line spectra will 
not be sufficient to fully determine the cloud map. However, once again MEM 
may be employed to locate the "smoothest" and "most symmetric" cloud maps 
that fits the data. Computers are now fast enough to support this type of 
detailed modelling and mapping of AGN emission regions. 

5.2. Mapping a Spherical Shell from Simulated Data 



Fig. 10 shows a simulation test designed to determine how well the geometry 
and physical conditions may be recoverable from emission-line reverberations. 
The adopted BLR model places clouds with density nn = 10 11 cm -3 and column 
density Nh = 10 23 cm -2 on a thin spherical shell of radius R/c = 12d. Synthetic 
light curves are computed showing reverberations in 7 ultra-violet emission lines, 
using CLOUDY to calculate the appropriate anisotropic and non-linear emissiv- 
ities. The synthetic light curves are sampled at 121 epochs spaced by 2 days, 
and noise is added to simulate observational errors. 

The MEM fit does not assume a spherical shell geometry, but rather it 
considers every possible cloud map *&(R, 9,nj{), and tries to find the simplest 
such map that fits the emission line light curves. The MEM fit adjusts 4147 pixels 
in the cloud map ^(R, 9, Tin), 143 points in the continuum light curve C(t), 7 
emission-line background fluxes L(X), and 1 continuum flux C. Note that the 
fit assumes the correct column density and distance. The fit to N = 968 data 
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points achieves x 1 /N = 1. The entropy steers each pixel in the map toward its 
nearest neighbors, thus encouraging smooth maps, and toward the pixel with 
the opposite sign of cos 6, to encourage front-back symmetry. 

On the left-hand side of Fig. 1C , we see that the recovered geometry displays 
the appearance of a hollow shell with the correct radius, and that the density- 
radius projection of the map has a peak at the correct density. The shell spreads 
in radius by a few light days, with lower densities at larger radii to maintain the 
same ionization parameter. 



On the right-hand side of Fig. 10, we see that the 8 light curves are well 



reproduced by the fit. The highs and lows are a bit more extreme in the data - a 
common characteristic of regularized fits. To the left of each emission- line light 
curve, three delay maps are shown corresponding to the maximum brightness 
(solid curve), minimum brightness (dashed), and the difference (dotted). All 
the lines exhibit an inward anisotropy except C in] . All the lines have positive 
linear responses except Mgn, which has a positive response on the near side of 
the shell and a negative response on the far side. 

This simulation test has used reverberation effects in emission-line light 
curves to map the geometry and density structure in the photo-ionized emission 
line zone. The 2-dimensional map clearly reveals the correct hollow shell geome- 
try, and moreover the correct density is also recovered. These promising results 
suggest that there are good prospects for probing the geometry and physical 
conditions in real AGNs. 



6. Two New Direct Routes to Hq 

Because the primary aim of this meeting is to discuss ways in which quasars 
may be useful as standard candles for cosmology, I will conclude with a brief 
discussion of two methods, both relatively new, which may allow the determina- 
tion of redshift-independent distances to active galaxies. Both methods employ 
echo mapping time delays to measure the physical size of something. Emission 
line time delays measure the sizes of regions of ionized gas surrounding the nu- 
cleus. Continuum time delays measure the sizes of reprocessing zones on the 
reverberating surface of an irradiated accretion disc. 

6.1. Hq from Sizes of Photo-Ionization Zones 

The ionizing photon luminosity Q can be estimated from an observed flux F as 

Q = AttD 2 SF , (31) 

where D is the distance and S is a "spectral shape" factor that converts the 
observed flux F to an ionizing photon flux 

SF =[ X0 M^. (3 2 ) 
Jo he 

If the ionizing radiation is not isotropic, the flux F in our direction differs from 
that seen by the photo-ionized gas clouds. That uncertainty may be absorbed 
into the factor S. 
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From the definition of the ionization parameter U, the radius of the photo- 
ionized zone is given by 



R= (^V>' = D (SFV*. (33) 

\4irUnftcJ \Uri}icJ 

Reverberation results allow us to measure R = tc from a time delay r. We 
therefore find that the distance is 

Since r is measured by reverberation, while U, nn, S, and F are found from 
the emission-line spectrum, we have a means of estimating the distance without 
using the redshift. For redshift z « 1, the Hubble constant is then 

Ho= D = t{u^-c) ■ (35) 

The above results apply only for the over-simplified case of photo-ionized 
gas characterized by single values of the ionization parameter U and density tih- 
This is intended only as an illustrative sketch of the method. In a more realistic 
model, the gas clouds have not single values but rather a cloud distribution 
function ^(R, 6, nn, Njj). We discussed this type of fitting in Section |5l. 

6.2. Hq from Reverberating Accretion Discs 

Steady-State Accretion Discs. The effective temperature on the surface of 
a steady-state accretion disc decreases with radius as 



/3GMM\ 1/4 
^ 8-KaR 3 ) 



(36) 



where M is the black hole mass and M is the accretion rate. If the disc sur- 
face radiates as a blackbody, the spectrum arising from the disc is obtained by 
summing blackbody spectra weighted by the projected areas of the annuli, 




D 2 

Here D is the distance, i is the inclination of the disc axis relative to the line of 
sight, and I = J °° x 5 / 3 /(e x - 1) ~ 1.932. 

The optical/ultra- violet spectra of AGNs have "Big Blue Bump" compo- 
nents that are attributed to thermal emission from accretion discs. Observed 
spectra are generally redder than the characteristic f u oc v 1 / 3 spectrum predicted 
by disc theory (e.g. Francis et al. 1991). The spectral signature of an accretion 
disc is therefore not very convincingly demonstrated. However, observed spec- 
tra are contaminated e.g. by starlight from the host galaxy (e.g. Welsh et al. 
1999). Taking difference spectra cancels out the contamination and measures 
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the variable component of the light. The difference spectra are usually bluer and 
can be in satisfactory agreement with the predicted spectrum for the change in 
brightness of an irradiated accretion disc (Collier et al. 1999). 
Echo Mapping Accretion Disc Temperature Profiles. In order to map 
the temperature profiles of accretion discs, we need some way to measure the 
temperature, and some way to measure the radius at which that temperature 
applies. We use a time delay to measure the radius, and a wavelength to measure 
the temperature. 

The disc surface is irradiated by the erratically variable source located near 
its centre, launching heating and cooling waves that propogate at the speed of 
light outward from the centre of the disc. This effectively reprocesses the hard 
X-ray and EUV photons that irradiate the disc into softer ultra-violet, optical, 
and infra-red photons. A distant observer will note that the reprocessed light 
arrives with a time delay 

R 

t = — (1 + sinicos60 . (38) 
c 

The heating wave requires a time R/c to reach radius R, and the reprocessing 
site at radius R and azimuth is at a distance from the observer that is larger 
by R sin i cos 9 relative to the centre of the disc. 

Note that for the annulus at radius R, the mean time delay, averaged around 
the annulus, is always R/c regardless of the inclination, but the range of time 
delays depends on the inclination. If the disc is face on (i = 0) then all azimuths 
have the same time delay. If the disc is edge on (i = 90°) then the far edge of 
the disc at 6 = has the maximum time delay r = 2R/c, and the near edge at 
9 = 180° has a time delay of zero. 

We use the blackbody spectrum to associate each wavelength with a tem- 
perature. The Planck spectrum for a blackbody temperature T, 

B »V) = -WW^T) ■ < 39 > 

peaks at the dimensionless frequency X = hc/XkT ~ 2.8. When irradiation 
increases T slightly, the change in the spectrum, proportional to dB u /dT, peaks 
near X ~ 3.8. 

To map the temperature profile of an irradiated accretion disc, we measure 
the time delay at different wavelengths. When we observe a change in the 
spectrum at wavelength A, we are measuring the reprocessed light from disc 
annuli where temperatures are 

T ~£x < 40 > 

with X ~ 4. Ultra-violet, optical, and near infra-red light curves give time 
delays for hot gas with T ~ 10 5 K, warm gas with T ~ 10 4 K, and cold gas with 
T ~ 10 3 K, respectively. If the temperature decreases with radius, the time delay 
should increase with wavelength. 

For the temperature profile of the steady-state disc, the time delay is 

21 



This r oc A 4//3 prediction has recently been verified in the case of the Seyfert 1 
galaxy NGC 7469 (Collier, et al. 1998). The observed r(A) therefore allows a 
measurement of MM, 

MM = (XX)-' r* . (42) 

Substituting for MM in the expression for f u , we find the distance, 

D = (^) 1/2 I 1/3 (^P) V2 (XX)- i/3 r . (43) 

Note that this distance is determined independently of the redshift of the AGN. 
Finally, for redshift z << 1, the Hubble constant is 

The above is of course only an outline description of the method. In practice 
one fits a model of reverberations in an irradiated disc to observed light curves. 
This new method is based on fairly straightforward physics - light travel time 
delays to measure radius, and blackbody spectra to associate a temperature 
with each wavelength. The results for NGC 7496 give Hq^Jcosi/^.I = 42 ± 
9 km s^Mpc" 1 (Collier, et al. 1999). 

The unknown disc inclination may not be too serious a problem because at 
high inclinations the dusty torus should obscure the broad emission line region. 
Seyfert 1 galaxies are therefore expected to have 0.7 iS Vcosi < 1. The v / cos7 
uncertainty may be reduced further by applying the method to a sample of 
Seyfert 1 galaxies with random inclinations in the above range. The use of 
blackbody spectra may be questioned, but can be checked by seeing if the same 
distance is found at different wavelengths. This test holds up for NGC 7469 
(Collier, et al. 1999). If the method can be shown to work for a larger sample 
of objects, it may serve to calibrate AGNs as standard candles for cosmology. 

Acknowledgments. Thanks to Brad Peterson and Shai Kaspi for provid- 
ing figures [H and ^ respectively. 
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